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Introduction 


e Many engine power-loss events reported since the 1990’s 


e Mason et al. hypothesized how power-loss events can result from 
ice crystals entering the engine core 


e Ingestion of ice into engine is studied at NASA PSL and elsewhere 
¢ Observed environmental conditions changed with cloud activation 
— Gas temperature change 
— Humidity change 
¢ Hypothesis: Thermal interaction between air and cloud 
¢ Model previously written to simulate NRC RATFac 
e Objective: Understand the air - cloud interactions in PSL tunnel 
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Model Formulation — General Description 


¢ Simulates PSL icing tunnel 


¢ Model couples air and cloud particle conservation eqs 
— Mass, energy fully coupled 
— Air is treated as ideal compressible gas 
— Isentropic equations used to solve pj; Vairs Tair 
e Full particle size distributions used 


e “air? = humid air = air + vapor 


Pp 
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Model Formulation - Assumptions 


¢ Air and particle flow are steady and one dimensional 

e Dry air and water vapor are ideal gases 

e Air (air + vapor) is well mixed 

¢ Tunnel is adiabatic and mass is conserved 

e Particle size distribution is characterized by a discrete set of diameters 

e Particles are evenly spaced 

¢ All particles are perfectly spherical 

e Particle aggregation and breakup through collision are negligible 

¢ Particles are injected in the direction of the flow and remain entrained 
¢ Temperature is uniform within the particle 

e Mixed phase particles are spatially homogeneous in water/ice content 

e Evaporation, condensation occur at the particle surface at particle temperature 
e The flow of particles and gas is a continuous stream 
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Model Formulation — PSL Description 


Tunnel Controllability 
e +0.3 kPa (.05 psia) 

* £05 °C (1.°F) 

¢ +1% RH 


8.84 m (29.0 ft) 


Inlet plane 


: —S Free Jet 
—— 

Spray Bars Exit Plane 
(Test Section) 
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Model Formulation — Experiment Configurations 
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Model Formulation — Differential Expressions 
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Model Formulation — Differential Expressions 


Gas Cloud Geometry 
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Model Formulation - Algorithm 


e Written in MATLAB version R2015b 
¢ Solves conservation differential equations using built-in ODE45 solver 


¢ Numerical relative and absolute convergence tolerance of 10° 
¢ Mass transferred between the gas and particle(s) balanced to 10°1° 


¢ Energy transferred between the gas and particle(s) balanced to 10% 
— Physical accuracy dependent on accuracy of property values (C,, Lyeat, etc.) 
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Sample Simulation 


Test Conditions 


To (K) 277.3 


RHo (%) 10.8 ‘TP#672| - 


Po (kPa ) 87.5 


——-Mixed Mass 


7 
: 


Ve (m/s) 85 


[JTunnel 


TWCrun(gim*)| 2.3 
Twater (K) 280.4 


MVD (um) | 42 


Tunnel Radius, m 
Particle Mass Fraction 


Twbo (K) 269.3 


Mixed Mass Ratio, g/kg 


0 5 10 


Tunnel Axial Distance, m 
Model Results 


Res. Time (S) | 1.26 
AToe (K) -2.7 


Mixed Phase 


ls 3 pm 


Awe (G/kg) | 1.2 


Me (melt ratio) | 0.20 = 40 pm 


&-330 wm 


TWC, (g/m*) | 1.0 —Total 


ATwhoe (K) 0.5 


Temperature, K 


Melt Ratio, Liquid/T otal 


5 } 5 10 
Tunnel Axial Distance, m Tunnel Axial Distance, m 
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Model/Experiment Comparison — TWC,,,, Sweep 


Ve (m/s) 
Po (kPa) 
To (°C) 
RHo (%) 
MVD (ym) 
Tater (°C) 


#672 


#671 
#670 y 
2.3 


0.78 1.4 
TWC, g/m? 


Melt Ratio, Liquid/Total 


AHumidity, g vap/kg dry air 


Takeaways: 
e Bulk vs Point 


© AWogaire ~ 30% 

© AT oaipe ~ 30% 

¢ ATwbo.= 
slight increase 


ATemperature, °C 
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Model/Experiment Comparison — RH Sweep 


Target Conditions ’ #667 @ Exp 
X 
‘ahs Sim Sim 
#669 
10 


Ve (m/s) 
Relative Humidity, % Relative ae % 


Po (kPa) 
To (°C) 


MVD (yum) 
Tater (°C) 


Melt Ratio, Liquid/Total 


AHumidity, g vap/kg dry air 


Wet Bulb Temps 


| 
(°C) | CC) , : 
Takeaways: ; yi 


¢ Twb, important 
e Small RH 
window 


ATemperature, °C 


10 35 
Relative Humidity, % Relative adhe % 
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Model/Experiment Comparison — 7,,,,., Sweep 


Target Conditions 
#667 #690 
@ Exp 
Bw Sim 
a 0 we 0 os 0 
43 82 


Ve (m/s) 
* Weker Temperature, °C Water ea perature, a 
Takeaway: 


Po (kPa) 
. @ Exp 
; wSim 


RAo (%) 


Melt Ratio, Liquid/Total 


TWCouik (g/ m?) 


AHumidity, g vap/kg dry air 


MVD (ym) 


ATemperature, °C 


To (°C) 
agreement o es 


43 82 
Water Temperature, °C Water Temperature, °C 
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Model/Experiment Comparison — Particle Size 
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Experimental Median Volumetric Diameter, pm 


Takeaway: Good MVD, agreement for MVD, = 15 um 
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Summary 


e Model written to understand Air - Cloud interactions in PSL 


¢ Model predicts to within 30% of measured changes in humidity 
and temperature 


¢ Model predicted satisfactorily for melt ratio 
— Some disagreement for elevated T,,.:., tests 


¢ Good agreement with particle size measurements 
¢ Twb, slight increase, important to determine cloud phase 


¢ Model guided development of test matrix for Fundamental 
Physics of ICI 2016 tests 
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